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Abstract. Non-holonomic mechanical systems can be described by a degenerate almost- 
Poisson structure @ (dropping the Jacobi identity) in the constrained space. If enough 
symmetries transversal to the constraints are present, the system reduces to a nondegen- 
erate almost-Poisson structure on a "compressed" space. Here we show, in the simplest 
non-holonomic systems, that in favorable circumnstances the compressed system is con- 
formally symplectic, although the "non-compressed" constrained system never admits a 
Jacobi structure (in the sense of Marie et al. 00). 



1. Introduction 

We adopt in this work the view of writing the equations of a non-holonomic systems in 
terms of an almost-Poisson bracket on a constrained manifold P C T*Q, introduced by 



van der Shaft and Mashke jlQfl . See also Cantrijn et al. [[3] for more recent developments. 

In this note we add a new twist to the simplest example of a non-holonomic system, 
the contact system in Q = R 3 (see eg. 0), namely: 

After performing the reduction by the transversal R 1 — symmetry, we observe that, for 
some metrics, the reduced ("compressed") almost-Poisson bivector admits a conformal 
symplectic structure, ie, a special Jacobi structure H0. In contrast, the "non-compressed" 
constrained system never admits a Jacobi structure. 

The examples point to the fact that, in quite favorable circunstances, the reduced 
system can be studied by symplectic techniques (for instance, when internal symmetries 
are present, integrability can be achieved by Marsden-Weinstein procedure, which holds 
in the conformally symplectic setting, see ||). The examples also show that, generally, 
non-holonomic systems are non- Jacobi systems (in the sense of Marie et al. 00). 

Hopefully our observations can help attracting more interest to investigations on the ge- 
ometrical properties of almost-Poisson bivectors which naturally describe non-holonomic 
dynamics on the constrained submanifold of the original cotangent bundle. 

2. The contact non-holonomic system 
Consider a non-holonomic systems in Q = R 3 having the constraint 
(1) z — xy = . 

The admissible sub-bundle E is the union of the horizontal spaces for the connection 1- 
form to = dz — xdy on the (trivial) bundle G = R 1 ■=— > Q = R 3 — > S = R 2 , with curvature 
duj = —dx A dy, where the R 1 action is the usual translation on the ^-fibers. 

For motivation, consider the javelin, a rod of mass m and moment of inertia I moving 
on a vertical plane (y, z) in such a way that it always remains tangent to its trajectory. If 
ip is the angle with the horizontal, then dz = tan(yj) dy, and if we introduce the change of 
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variables x = tamp, then T = ~ [m(y 2 + z 2 ) + lip 2 } = ~ [m(y 2 + z 2 ) + I n^ay ] • Choose 
units so that m = I = 1. Under the assumption of small angles, (p ~ x, we get 

(2) L = r = i(x 2 + 2 / 2 + i 2 ). 

According to our previous work (0) system is "z-Caplygin" and thus can be reduced 
to TS = {(x, y, x, y)} with lagrangian 

(3) L= l -{x 2 + {l+x 2 )y 2 ) . 

with an external gyroscopic force added to it. Actually it can be also interesting to 
regard it as a y-Caplygin system dy = cot(^) dz = (l/x)dz , since we may want to add 
the gravitational potential V = gz. 

3. Almost-Poisson brackets via moving frames 

The giroscopic force can be concealed in an almost Poisson bracket in the constrained 
manifold P C T*Q via the dynamical equation 

x = {x,H} P 

where H is the hamiltonian, and P C T*Q is the Legendre transform of the constraint 
subbundle E C TQ, as defined below. 

Let ej = eijd/dqL be a moving frame on Q such that the first m vectors (labelled by 
lowercase latin indices) generate E. By a direct calculation ([p!0[, equation(19)) van der 
Shaft and Maschke verified that the brackets are given by 

{Qi,Qj} = °) UiiPj} = e L j(q) = dqi ■ ej (ej = e LJ d/dq L ) 

(4) {Pi,Pj} = ~Pq ■ [ej, ej] = R u 

where p q is evaluated on P. For a geometric interpretation and simple derivation of these 
formulas using the moving frame method see our paper ]/]]. Here's an outline: 

Choose an adapted moving frame to the constrained distribution E C TQ that is, 
consider a complete set of vector fields e« , e a where e,(g) G E q . The greek labels a denote 
e a ^ E g . Denote the dual 1-forms by ej , e a . We shall denote by uppercase latin indices 
ej , Si the complete dual set in TQ and T*Q. The canonical 1-form on T*Q writes as 

(5) pdq = ujej 
where ui define new coordinates on each T*Q. Now, 

d(pdq) = duj A£; + ujdei . 
The last term vanishes on vertical vectors. Moreover 

u/de/(ej,e^) = U/e J (e I (e i t)) - u/e^(e 7 (ej)) - u 7 ej[ej, e K ) = -ti/e/fej, e K ] 
which gives us the symplectic matrix on T*Q as 

(6) [0] 

moving frame I j r\ 

and its inverse, the Poisson matrix on T*Q as [A] moving coframc = [^] m ovin g coframc > where 

(7) = -w/e/[ej, ejf] = -p 9 • [e J? e K ] . 
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The moving coframe for [A] movingco f rame is ej,duj. The moving frame for is its dual. 
Caveat: this basis contains the vertical vectors djduj and lifted e} = tj + vertical (see 
0). From this we obtain the matrix of almost-Poisson brackets in P C T*Q by directly 
substituting the u a by the Ui (and q's) via the defining equation for P : 

(8) dH/du a = 

which is a consistency requirement for any constrained (holonomic or non-holonomic) 
dynamical system on Q , see 0. 

4. The contact almost-Poisson structure 
In the case of the contact system, let us begin by taking as moving co-frame the set 

(9) ei = dx , e 2 = dy , €3 = to = dz — xdy 
that is dual to 

(10) ei = d/dx , e 2 = d/dy + xd/dz , e 3 = d/dz . 
Note that 

(11) [ei, e 2 ] = e 3 , [e u e 3 ] = [e 2 , e 3 ] = , 

so that we have the 3-dimensional Heisenberg algebra. The identity @ yields 

(12) Ui = Pa , U 2 = Pj, + Xp x , U 3 =p z , 

hence the Hamiltonian is 

(13) if = ^(pI+pI + pI) = 1 -{ U \ + (n 2 - xm 3 ) 2 + u^) . 

Now, the Legendre transform of the admissible sub-bundle E C TQ is given by (|) as 
dH/dus = so that the constrained manifold P is 

XU2 

(14) P = {(x,y,z,U!,u 2 ,u 3 ) \ u 3 = — — r } = {(a:, y, |p z = xp y } . 

1 + x z 

As moving basis for T*P we take <ix, dy, dz — xdy, du\, du 2 and the corresponding 5x5 
matrix of Poisson Brackets is given by 







f o 2x2 


2 xl 


^2x2 


(15) 


[A] moving 1 


0lx2 


Oixl 


0lx2 






V —12x2 


2x i 


-^2x2 



where R 2x2 is the antisymmetric matrix with 

(16) Ru = -p q [ei, e 2 ] = -u 3 = - ■ 

Here, the vanishing of the middle row and column means of course that 

(17) {dz — xdy,H} = z — xy = 
which gives the differential equation for z. 
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5. The compressed system 
Deleting the middle column and row of fll5D , we obtain a non-degenerate matrix 

TaT — r a i — I ^ 2x2 ^ 2x2 

I/ 1 - J l 1 *-J compressed I r r> 

\ — J 2x2 -"-2x2 

which characterizes the almost-Poisson structure in P = P C om P ressed = {(x, y, ui, u 2 )} with 



(18) 



(19) 



H = H, 



compressed 



2 lU 



11.-, 



ar 



as the reduced Hamiltonian. It follows that the reduced equations of motion are 



or 



x 



ui , y 



X 

y 

U 2 
U 2 



02x2 
~hx2 



1 2x2 
R2x2 



H_x 

Hy 

H_U! 

H 



u 2 



1 + X 2 



Ui 



-Hr + R\ 2 H 



U 2 



—H y — Ri 2 H Ul 



A "miraculous" cancelation takes place in the iii equation, and the system is actually 

U2 . „ . x 



(20) 



x = m,y 



1 , 2 ' Ul 
1 + ar 



, u 2 



;UiU 2 



1 + x 2 

This cancelation has a reason. The Lagrangian is invariant under the 1-parameter group 
(x, y, z) — > (x + e,y, z) and the generator d/dx is an admissible vectorfield. By the non- 
holonomic Noether theorem [0, the momentum p x = U\ = x is conserved. The 2 degrees 
of freedom compressed system separates and can be integrated by quadratures: 



ln(w 2 ) 



x + at 



U2 



1 + (x + at) 2 
Ay/1 + (x + at) 2 



dt 



- ln(l + (x Q + at) 2 ) + const. Thus, 



y = y + 







1 + (x + at) 2 

and we reconstruct the z-fiber dynamics via the constraint equation: 

'* ^ Ay/l + (x + at) 2 ^ 
(x + at) — ; ; — — at . 



z + 



1 + (x + aty 



The almost-symplectic form in the compressed space T*S, S 



R 



{x, y}, is the 



2-form Q = du\ Adx + du 2 Ady + R\ 2 dx A dy , R\ 2 = — and so dfl = —jfp: du 2 dxdy . 

Let us investigate if there is a function / : T*S — > R such that fCl is closed. It looks 
simpler to try the ansatz f = f(x) only. The condition df AQ + fdfl = leads to 



(21) 



if(x) + f 



x 



1 + x' 



:}dxdu 2 dy = , so that / 



.4 



The significance of this observation is that the compressed system (|20D is Hamiltonian in 
the time scale s such that dt/ds = a/1 + x 2 , with the same Hamiltonian H = H compresse d 
and bona-fide symplectic form fQ. Now, the corresponding Poisson structure in T*S is 
j [A] compressed (refer to (p^8|)) and one is tempted to guess that the conformally changed 
bivector in P given by y[A] moving (refer to fll^D) satisfies Jacobi. But as we shall see below, 
this is not the case / 
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6. Non- Jacobi for the constrained almost-Poisson 

Let us now take a closer look on the algebraic properties of the almost-Poisson structure 
on the constrained space P. In order to compute the Schouten-Nijenhuis bracket of the 
almost-Poisson bivector A with itself, we first rewrite A given by (]T5|) as a matrix in a 
coordinate basis. Choosing {x, y, z, u\, 1*2} as coordinates, we write A^ Xt y jZjUl)Ua y as 



02x2 


02x1 


h 


0lx2 


Olxl 




—12x2 


-T T 

-^2x1 


-^2x2 



"2 



where, as before, i?2x2 is the antisymmetric matrix with R12 = —xv/(l + x 2 ) and now 

L lx2 = ( x ) 

in such a way that we can compute the self S-N bracket [A, A] using: 

(22) [A, A] IJK = A LK d L A IJ + A LI d L A JK + A LJ d L A KI 
where the summation convention is subtended. We then get 

(23) [A, A] = (j^J (xd y - d z ) A d Ul A d„ 

One can check explicitly that there is no vector field E satisfying the first of the equations 
for the existence of an associated Jacobi structure |3| (the second one is given by (|4"0D) : 

(24) [A, A] = 2E A A . 

Thus no Jacobi structure, in the sense of Marie, exists for this bi- vector (or any conformal 
one as well). Notice that we can rewrite [A, A] = [jt^s) Ad Ul Ad U2 , where e 1 - = xd y — d z 
is a vector orthogonal to e\ and e 2 and hence to the distribution. Remark that A// , 
where / is given by (|2ID, is really "^-almost" Poisson. A calculation in the same lines 
shows that the only nonvanishing entries of [A//, A//] are the permutations of (3,4,5), 
namely —2(1 + x 2 )d z A d Ul A d U2 . ("almost" is really almost)! 

In order to better appreciate the Non- Jacobi result for the almost-Poisson structure, let 
us now choose a different moving frame and co-frame adapted to the contact distribution. 
Specifically, we choose an orthonormal set with respect to the euclidean metric: 

(ok) P -a p - dy + xdz P - dz ~ xdy 

(25) ei - d x , e 2 - , e 3 - 

Vl + i v 1 + x 2 

whose dual set is "itself" : 

dy + xdz dz — xdy uj 

(26) ei = dx , e 2 = , , e 3 



VTTx 1 VYTx 1 VTTx 1 ' 

This frame does not respect the natural z-fibration, but it's still true that cj(ei) = co>(e 2 ) = 
0. However, the Lie algebra is now modified : [ei, e 2 \ = (j^s) e 3 5 [ e 2, = , [e^, e\] = 
(1+^2 ) e 2 > but the hamiltonian H corresponding to the euclidean metric is again euclidean: 
H = (uf + u\ + w|)/2 , so that the condition dH/du?, = yields the simpler equation 

(27) u 3 = 

for the definition of P (generally, we have u a = for orthonormal frames). 
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From the theory of moving frames M, the Poisson bi- vector on T*Q can be written as 

(28) A = e} A d Ul + Rud Ul A d Uj 

where Ru is given by ([|), d Ul is the "vertical dual" to dui and e} is the correct lift of 
the base vector ej to T(TQ) M. Alternatively, we can rewrite this full bi- vector using 
Darboux quasi-coordinates as A = ej A d Ul . For an orthonormal adapted moving frame, 
from (5) we get that the almost-Poisson bivector on the constrained space P is: 

(29) A = ei Ad Ui + Rijd Ui A d Uj 
where ~ means evaluating at P, in this case: 

(30) R u = -«iej[e/,ej] . 

For the contact system with orthonormal moving frame, we have simply 

(31) A = e x A d ui + e 2 A d U2 

where d Ui is the "vertical dual" to ej : d Ul = d Px , d U2 = 9p ^jz^ z > d U3 = 9v ^^ v ■ The 
self S-N bracket of A can be easily computed using the decomposition formulas for the 
S-N bracket of bivectors. We get : 

[e 1 Ad ul +e 2 Ad U2 ,e 1 Ad Ul +e 2 Ad U2 } = [e x Ad Ul , exAd Ul ] + [e 2 Ad U2 , e 2 A<9 U2 ]+2[e 1 A<9„ 1 , e 2 Ad U2 ] 
It's easy to see that the first two terms on the right vanish. The third one decomposes as 
[ex Ad Ul , e 2 Ad U2 ] = [ex, e 2 ]A d Ul A d U2 - e x A [e 2 , d Ul ] A d U2 - e 2 A [e x , d U2 ] A d Ul + e x A e 2 A [d Ul , d U2 ] 
Since the only nontrivial Lie bracket is [ex, e 2 ] = [jt^i) &3 we are left with 

(32) [A, A] = (r^j e 3 A d Ul A d U2 , 

which is equivalent to the result previously obtained (^). It's now even simpler to verify 
that no Jacobi structure [f||| is possible for A, just compare ([24]), (|31|) and ([32]). 

Let us consider a still simpler, even more symmetrical case. Namely, on the con- 
tact system, we consider a metric which is invariant under the full Heisenberg group. 
Again, taking (|10D as moving frame we have the Heisenberg algebra ([lTJ) . The simplest 
Heisenberg- invariant metric is thus given by the kinetic energy T = ^(vf +v\ + vf) , where 
Vi is given by the identification Vxtx + v 2 e 2 + v 3 e 3 = xd x + yd y + zd z , which gives: 

(33) T = ^(x 2 + (1 + x 2 )f + z 2 - 2xyz) . 

In other words, for such metric 2T, the Heisenberg moving frame (|10|) is orthonormal. In 
terms of the hamiltonian, we have that the co-frame @ is also orthonormal, and thus 

(34) H = ^{u{ + u\ + u\) = -{pi + v\ + (1 + x 2 )pI + 2x PyPz ) 

where the relation between the u's and the p's is the one given earlier in (0), section 3. 
Once again, the almost-Poisson bivector is given by 

(35) A = ex A d Ul + e 2 A d U2 , and thus [A, A] = 2e 3 A d Ul A d U2 . 
Therefore, no Jacobi structure ||@ exists on P for this Schaft-Maschke [|It| bivector. 
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7. The compressed almost-Poisson structure 
is not conformally symplectic in general 

We've seen earlier in section 5 that the compressed system for the contact nonholo- 
nomic system with euclidean metric is conformally symplectic. However such conformal 
symplectic structure (or any Jacobi structure) on the compressed system does not exist 
for generic metrics. To see this, consider a ^-invariant hamiltonian of the general form 

(36) H = jijUiUj + V(x, y) , j tj = j tj (x, y) = j^x, y) . 

The condition dH/du^ = yields = — (713W1 + 723W2V733 = —R12 > so that the 
compressed almost-Poisson bivector can be written as 



(37) 



A = d x a d ui + d y a d U2 + ( lnUl + l23U2 ) d Ul A d U2 . 

V 733 / 



Since {x, y,Ui,u 2 } is a coordinate system for T*S, we can apply formula (|22|) directly to 
obtain [ A, A ] = -2 (^j d x A d Ul A d U2 — 2 (^j d y A d Ul A d U2 . On the other hand, a 
general vector field on T*S has the form : 

(38) E = ad x + Pd y + fxd Ul + vd U2 
so that the condition [ A, A ] = 2E A A is satisfied for 

(39) v = -713/733 , H = 723/733 • 

For E to be the vector field of a Jacobi structure || , it is also necessary that 

(40) [E,A] = C E (A)=0 . 

A simple computation gives [E, A] = — (j^ {j^j + jfj (^JJ ®ui A d U2 and thus 

9 ( 713 \ d ( 723 \ 
41 dx V733J + dy V733J ° 

is a necessary condition for the existence of a Jacobi structure on the compressed system 
in T*S. Of course, the hamiltonians (0) and (|34D , obtained from the euclidean and 
the Heisenberg metrics, (§) and ( |3~3D respectively, satisfy the above condition. More 
generally, in order to verify this condition, one can substitute for the original metric 
elements gij(x\, 22) in a basis x\ = x, 22 = y, 23 = z, with the relations 

^ 713 _ 9l2g23 - 913922 ~ Xx {9X2933 ~ 913923) 



(43) 



733 9u922 ~ g\ 2 - 2x1(5-115-23 - #12513) + x\{9n933 - 9I3) ' 

723 = 011023 ~ 012013 ~ ^1 (511533 ~ gfgj 

733 011022 - 0? 2 - 2Zi (0H023 - 012013) + A. (011033 ~ 0? 3 ) ' 



Clearly, (|3"8D, ( |3T3| ) and (O) mean that any possible Jacobi structure on this compressed 
system comes from a conformal symplectic structure. 
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8. Conclusions 

We've seen that the failure of the Jacobi identity for the Schaft-Maschke almost 



Poisson structure on the constrained submanifold P C T*Q of a nonholonomic system 
cannot be cirvumvented by the introduction of an associated Jacobi structure ||] || , even 
in the simplest cases. It remains to be confirmed whether this "no-go" result is generic. 

These examples also suggest that pehaps the almost-Poisson structure can be most 
useful when there is a principal bundle structure G — > Q — > S, and all data are equivariant 
with respect to the group G. In the contact case, z-invariant Lagrangians. In this case, 
the degenerate almost-Poisson bracket in P projects over a non-degenerate almost-Poisson 
bracket in T*S, which in favorable cases (not all) is conformally symplectic, or Jacobi. It 
remains to be found a geometrical interpretation for these cases. 

More generally, it also remains to be determined whether or when such a conformally 
symplectic reduction can be applied for less simple nonholonomic systems. Work in this 
direction is under way and shall be reported elsewhere. 



[1 
[2 
[3 

[4; 

[5 
[6 
[7 
[8 
[9 

[10 



References 

V.I. Arnold, V.V. Kozlov and A.I. Neishtadt, Mathematical aspects of classical and celestial 
mechanics, Dynamical systems. Ill , Encyclopaedia Math. Sci., 3, Springer, Berlin, 1993. 
L. Bates AND R. Cushman, What is a completely integrable nonholonomic dynamical system? 
Rep. Math. Phys., 44: 1/2 (1999), 29-35. 

F. Cantrijn, M. de Leon and D. Martin DE Diego, On almost-Poisson structures in non- 
holomic mechanics, Nonlinearity 12 (1999), 721-737. 

P. Dazord, A. Licherowicz and CM. Marle, Structure locale des varietes de Jacobi, J. Math, 
pures et appl. 70 (1991) 101-152. 

S. Haller and T. Rybicki, Reduction for locally conformal symplectic manifolds, J. Geom. Physics 
37:3 (2001), 262-271. 

J. Koiller Reduction of Some Classical Non-Holonomic Systems with Symmetry, Arch. Rational 
Mech. Anal. 118 (1992) 113-148. 

J. Koiller , P. de M. Rios and K.M. Ehlers, Moving Frames for Cotangent Bundles, to 



appear in Rep. Math. Phys. Eprint: math-ph/0203014 



M. DE Leon AND D. Marti'n de Diego On the geometry of non-holonomic Lagrangian systems, 
J. Math. Phys., 37 (1996), 3389-3414. 

CM. Marle, On Jacobi Manifolds and Jacobi bundles. In "Symplectic geometry, groupoids and 
integrable systems" (P. Dazord and A. Weinstein, eds.), MSRI Publications 20 Springer (1991) 
227-246. 

A.J. van der Schaft and B.M. Maschke, On the hamiltonian formulation of nonholonomic 
mechanical systems, Rep. Math. Phys. 34:2, (1994) 225-233. 



Laboratorio Nacional de CoMPUTAgAO Cientifica, Av. Getulio Vargas 333, Petropolis, 
RJ 25651-070, Brazil. Present address (Rios): Dept. of Mathematics, Univ. of California, Berkeley. 
E-mail address: prios@math.berkeley.edu , jair@lncc.br 



